Nucleation of plasmonic resonance nanoparticles 
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We predict the electromagnetic field driven nucleation of nanoparticles that provide plasmonic 
oscillations in resonance with the field frequency. The oscillations assume a phase that maximizes 
the particle polarization and energy gain due to nucleation. We derive closed form expressions for 
the corresponding nucleation barrier and particle shape vs. field frequency and strength, metal 
plasma frequency, conductivity, and the host dielectric permittivity. We show that the plasmonic 
polarization allows for nucleation of particles that would not be stable in zero field. 



PACS numbers: 

Accelerated nucleation in response to laser or dc elec- 
tric fields has been observed in a number of systems*^— 
The phenomenon is typically attributed to the lower- 
ing of the nucleation barrier by the field-induced po- 
larization of the new phase particle. Similar observa- 
tions have been reported for nanoparticles of various 
shapes.- - ^ Existing theoretical models^r— consider an 
essentially static field. However, for ac fields, especially 
near plasmonic resonances, frequency dependent effects 
can strongly modify the nucleus polarization and expo- 
nentially change the nucleation rate. 

This Letter introduces the concept of plasmonic medi- 
ated nucleation. The underlying idea is that the nucle- 
ated particle provides a narrow resonance and phase of 
plasmonic oscillations that maximize its polarizability in 
the external electromagnetic field; that makes such par- 
ticles energetically more favorable and lowers their nu- 
cleation barriers. As a result, an ac field will favor the 
nucleation of a particular geometry, such that the plas- 
monic frequency of the particle is in resonance with the 
field frequency. Furthermore, we will show that the plas- 
monic related energy gain can be significant enough to 
change the phase equilibrium, thereby triggering nucle- 
ation of particles that would not form in zero field. 

Our consideration is based on the classical nucleation 
theory (CNT)^i~— where the system free energy is ap- 
proximated by the sum of the bulk and surface contribu- 
tions of the new phase particle. Adding an electric field 
dependent term, Fe, the free energy takes the form, 



F = F E + fiV + a A. 



(1) 



Here, /x is the difference in chemical potential (per vol- 
ume) due to nucleation, and a is the surface tension, 
V and A are the particle volume and area respectively. 
The case of fj, < corresponds to a metastable system in 
which nucleation is naturally expected without external 
field; fj, > describes the case where metal particles are 
energetically unfavorable in zero field, yet, as shown be- 
low, they can appear in a sufficient electromagnetic field. 

For a static fielder— Fe represents the polarization 
energy gain similar to that of an induced dipole p = aE 
in an electric field E. Given a particle of polarizability 



a in a dielectric material with permittivity e, the polar- 
ization induced gain in a static field can be represented 
as^i 



-eaE z 



(2) 



Although not immediately obvious, the factor e makes 
Eq. ^ different from the energy of a dipole in an exter- 
nal field. It reflects the contributions from all charges in 
the system, including those responsible for the field; this 
factor has been confirmed by several authors^— 

In order to introduce some useful notation, we briefly 
review the case of a spherical metallic nucleus of radius R 
in a static field. Using a = R 3 , V = 4tt j R 3 /3, A = 4ttR 2 , 
and assuming /i < 0, the nucleation barrier W and radius 
R are determined by the maximum of the free energy in 
Eq. ©, 

W = W (l + £/2r 2 when R = R {1 + £/2) -1 . ( 3 ) 

The corresponding zero-field CNT quantities and the di- 
mensionless field strength parameter £ are, 
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(4) 



Their ballpark values are Wo ~ 1 eV, Ro ~ 1 nm, and 
( C 1, say £ ~ 10~ 5 for a moderate field of E = 30 
kV/cm and e = 1. 

Since we are considering electromagnetic fields of fre- 
quency ui 3> uj at , where u at ~ 10 13 s _1 is the characteris- 
tic frequency of atomic vibrations, the induced polariza- 
tion will be predominantly of electronic origin. The cor- 
responding part of the electronic energy, which is propor- 
tional to — p • E, will oscillate in time with frequency to. 
According to the standard procedure of adiabatic (Born- 
Oppenheimer) approximation in solids, the time average 
of that oscillating energy can be treated as a contribution 
to the potential energy of the atomic subsystem. Taking 
into account the known recipe for time averages^ the 
latter contribution to the free energy of an atomic sub- 
system can be presented in the form, 



Fe = 



(5) 



2 



where E is now understood as the field amplitude and 
5J(a) represents the real part of the polarizability. 

Prior to the rigorous derivation, we show next that 
all the major results of this work can be obtained qual- 
itatively [Eqs. (foTITUl)]. We start with noting the high 
polarizability, 



a ~ {H/R) 2 V > V, 



(6) 



of needle-shaped particles, e.g., a prolate spheroid or a 
cylinder of height 2H and radius R <C H aligned to the 
field (Fig. [1]). Indeed, the field- induced charges, ±q, 
induced at the opposite poles are estimated from the 
balance of forces, q 2 / H 2 = qE, which gives the dipole 
moment p ~ Hq ~ EH 3 ~ V(H/R) 2 E = aE. 




H 



FIG. 1: Illustration of a plasmonic resonance nanoparticle. A 
metallic prolate spheroid of semi-major axis H (aligned with 
the electric field E) and semi-minor axis R is formed with its 
aspect ratio tuned to the frequency, u>, of the ac field such 
that H/R w Ldp/uj, where ui p is the plasma frequency, e and 
e p are the permittivity of the dielectric medium and particle, 
respectively. 

The plasmonic resonance in a needle-shaped metal 
particle can be qualitatively explained by considering 
small parallel displacements x <C H of the free electrons 
along the cylinder axis. This will result in small charges 
q ~ ±irR 2 xNe on the two poles of the cylinder where 
N and e are the electron concentration and charge, re- 
spectively. The corresponding forces —qe/H 2 on individ- 
ual electrons can be interpreted as the restoring forces 
—rmo 2 x, which yields the resonant frequency, 



u r ~ co p (R/H) <C Up, 



(7) 



where m is the electron mass and we have used the stan- 
dard definition of the plasma frequency, 



uj p = ^jAirNe 2 Jr 



(8) 



oo r is known as the plasmonic resonance frequency which 
describes collective oscillations of quasi-free electrons; it 



has been experimentally observed in light scattering by 
nanoparticlesJ^ 

When uo r is in resonance with u), the plasmonic vi- 
brational energy is accumulated in the particle to the 
extent allowed by its quality factor (Q-factor), estimated 
as Q ~ lot 3> 1 with r being the electron relaxation time. 
Because that energy is proportional to the square of the 
internal electric field, Ef nt , one can say that the latter 
increases the polarizability effectively by the factor uj r r. 

Combining the latter with Eq. ^ and (J7J gives the 
resonant (to u r ) plasmonic enhanced polarizability, 



V(0J p /0j) 2 0JT > V. 



(9) 



This gigantic increase in polarizability takes us to the 
major prediction of this work: an ac field of frequency ui 
can drive the nucleation of needle-shaped particles with 
resonant aspect ratio H/R ~ uj p /uj 1. 

We note that the first amplification factor (w p /w) 2 ~ 
(H/R) 2 3> 1 reflects the needle-shaped geometry of the 
particle and remains as such in a static field^ while the 
Q-factor enhancement lot 3> 1 is specific to ac fields. 
Taking into account the amplification ratio in Eq. ([S]), an 
approximate result for the ac resonant nucleation barrier 
can be guessed from the known static result^ or even 
from Eq. (O, with £ — > (lo p /lo) 2 (lot)£, > 1, which yields, 

W ~ (lo/lo 3 p t 2 )(W /E 2 R 3 ) 2 W . (10) 

Our rigorous analysis begins with the polarizability of 
a spheroid)^ 



V 



4-7T e + n(e p — e) 



(11) 



Here, e p is the dielectric permittivity of a metal particle 
and e is that of the medium, n is the depolarizing factor, 
which we take as that of a strongly anisotropic prolate 
spheroid, 

n k (R/H) 2 [\n(2H/R) - 1] = {R/H) 2 A < 1. (12) 
Also, we use the dielectric permittivity of a metal, 



LO^ LO°T 



(13) 



e is assumed to be a real number, and r 1 -C oo <§C uo v . 



Using Eqs. ([TT1-H3]l with Eq. © yields, 
V (n — nj) + bn 



R(a) 



47r (n — n^) 2 + bn 2 ' 



where b = \/(lot) 2 and, 



+ IV 2 



< 1. 



(14) 



(15) 



The polarizability, %i(a), has a sharp minimum when, 

n k n u - Vbn^, (16) 



3 



which reflects the presence of the plasmonic resonance. 
Since n ~ (R/H) 2 from Eq. (TT2"1) . the expression for n u 
is fully consistent with the qualitative result in Eq. (0. 
The relationship between nanoparticle aspect ratio and 
the field frequency is illustrated in Fig. [2] 
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FIG. 3: Normalized free energy of a prolate spheroidal metal- 
lic particle vs. field frequency / = uj/2iy, with E — 3 x 10 5 
V/cm. The sharp resonance of width 1/r determines the as- 
pect ratio H/ii ~ uj v /uj ~ 10. Parameters values are typical 
for nucleation in solids: Rq — 1 nm and Wo = 1 eV, with 



= 10 ib rad/s, r = 10" 



and e = 1. 



in metastable systems. For the case of stable systems, 
the field can induce nucleation when. 



FIG. 2: Aspect ratio of the metallic nucleus as a function 
of frequency. The resonance frequency for a sphere is ui r = 
u p /V3, while for a prolate spheroid the resonance frequency 
depends on the aspect ratio, as shown. 

Given the sharpness of the resonance, all other n- 
dependent quantities, in particular the prolate spheroid 
volume and area, 



V = 



3A : 



A = ir 2 H\ ^ 



can be evaluated at n = n u . That yields, 



!»(«)], 



H A 1 



(17) 



(18) 



fully consistent with our earlier estimate in Eq. (|9]). 

Normalizing the free energy in Eq. ([1]) with respect to 
the classical barrier, it takes the form, 



E > E r 



1 



OJT 0J V 



iQWo 



(20) 



Assuming the above ballpark parameter values yields the 
characteristic field E c ~ 10 8 V/cm. However, the other 
multipliers can make the right hand side in the inequality 
of Eq. (|20|) much lower, down to say 10-100 kV/cm, 
corresponding to a moderate laser power density of P ~ 
10 - 100 mW/pm 2 . 

From this point on we concentrate on the case of field 
dominated nucleation, which satisfies the condition of Eq. 
(|20p and implies neglecting the second term in the square 
brackets of Eq. (|19[) . In minimizing F with respect to H 
we keep 



A«ln< 2uJp 



(21) 



constant. The resulting nucleus semi-major axis and nu- 
cleation barrier are, respectively, 



F 
Wo 



Rl 



2e 
1 



E 2 eRl 



24 W, 



3tt U u H 2 



(19) 



The minus sign in the second term corresponds to the 
case when the original phase is metastable and, thus, nu- 
cleation can occur in zero field. The plus sign describes 
the case of a stable original phase, i.e., nucleation is im- 
possible in zero field. The effect of the resonance on the 
free energy is illustrated in Fig. [3] 

It follows from Eq. (IT9l that the plasmonic reso- 
nance mechanism can significantly accelerate nucleation 
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i 



-Ra 



§) n-„ 



(22) 



(23) 



The above theory applies most easily to a system close 
to some phase transition at temperature T c . That enables 
one to estimate the chemical potential as [i = ^o(l ~ 
T/T c ). Correspondingly, the classical nucleation radius 
and barrier become, 



Rq — floo(l — T/T c 



W = W m (l-T/T c 



(24) 
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and E c = E cQ (l - T/T.) 1 / 2 , where R 00 , W 00 , and E c0 
are obtained from their definitions in Eqs. (U|) and ([20)1 
with [j, = [1q; numerically, they are of atomic scale. The 
convenience of this model is that it allows macroscopi- 
cally large Rq, consistent with CNT; also, it corresponds 
to lower a critical field, E c , making the plasmonic nucle- 
ation effect easier to observe. 

Using the above scaling, the results for particle nu- 
cleation length and barrier turn out to be temperature 
independent: they retain their form of Eqs. ([22]) and (|23|) 
with the trivial substitutions, 

Ro -> R 00 , Wo -> Woo, E c -> E cQ . 

This is in striking difference with CNT, which predicts a 
diverging nucleation radius Ro and barrier Wo towards 
the phase transition temperature in Eq. (|24[) . 

Numerically, there is a broad window wherein the 
above predicted plasmonic resonance nucleation takes 
place. Indeed, the quality factor at the plasma frequency 
can be represented^ as uj p t = 160 / (y/ Na^p) ~ 10 3 , 
where as is the Bohr radius. The factor 1/ ^Na? B is typ- 
ically in the range of 5-10, and the resistivity, p, (in units 
of p£l cm) is typically smaller than unity. However, this 
estimate may become less optimistic if r is additionally 
reduced by particle surface scattering when the particle 
dimension is very smallJ^ 

We now briefly discuss possible implications of our the- 
ory. A system near phase transition under a moderate 
intensity laser beam would be most suitable for experi- 
mental verification. The predicted high aspect ratio par- 
ticles are optically accessible and can be identified via 
the unique features of prolate spheroids in light scatter- 
ing and absorption . 15 ! 18 Another unique feature is that 
the laser beams that create the nanoparticles can simul- 
taneously play the role of optical tweezers. That can be 
exploited to observe and control the positions of nucle- 



ated metal particles that would be unstable in zero field. 

There have been several observations that are possibly 
consistent with the above predicted mechanism, a few 
of note include: the formation of metallic nanorods un- 
der laser beams, which did not nucleate in zero field; 5 
many orders of magnitude acceleration of metal nano- 
particle formation under laser irradiation;— laser induced 
nucleation of silver nanoparticles in glasses^ and a ver- 
satile, comprehensive study of laser induced nano-wire 
nucleation^ In all those cases, different mechanisms may 
have been responsible for the observations (except, per- 
haps, for the case 5 - that specifically identifies the fre- 
quency dependent effect and the possible role of plas- 
monic excitations). Further experimental verification is 
required to attribute the observations to the phenomenon 
predicted here. Verifiable features include frequency de- 
pendent nanoparticle nucleation and shape, nucleation 
rates that are exponentially accelerated with laser power, 
and Arrhenius temperature dependence of the nucleation 
rate. It should be noted that, because the laser beam ac- 
tion continues for a certain time, the nucleated particles 
can evolve further, undergoing growth, which can affect 
their expected shapes. 

Given the abundant current literature on the topic, 
we can safely omit the discussion of important practical 
applications of nanoparticle plasmonic resonances. How- 
ever, we emphasize that because of its resonant nature, 
the mechanism introduced in this work can be useful in 
better controlling nanoparticle sizes, shapes, and distri- 
butions. 

In conclusion, we have predicted a phenomenon of plas- 
monic driven nucleation of metallic nanoparticles. There 
are many important practical implications. More work is 
called upon to describe the growth stage of such resonant 
particles and relate them to the experimental observa- 
tions. 
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